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Abstract 

We extend the Gyulassy-Levai-Vitev reaction operator approach to multiple elastic scattering 
of fast partons traversing dense nuclear matter to take into account the leading power corrections 
due to the medium recoil and to derive the change in the partons' longitudinal momentum. We 
employ a boost-invariant formalism to generalize previous treatments of the problem, which were 
specific to the target rest frame. We apply the transverse momentum diffusion results in a simple 
analytic model to evaluate the broadening of the back-to-back di-hadron correlation function in 
d + Au reactions. 
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I. INTRODUCTION 



New experiments at the Relativistic Heavy Ion Collider (RHIC) will provide high statistics 
measurements of photons, leptons and hadrons at large transverse momentum (jpx) in high 
energy nucleus-nucleus collisions [lj. Particle production from a single hard scattering with 
momentum exchange much larger than 1/fm should be localized in space-time. It is multiple 
parton scattering before or after the hard collision that is sensitive to the properties of the 
nuclear matter fl Q. B y C „,„pa ri „ g the high-p, observable* fl fl Q in p + p, d + A 
and A + A reactions, we are able to study the strong interaction dynamics of QCD in the 
vacuum, cold nuclear matter and hot dense medium of quarks and gluons, respectively 

In this Letter, we extend the Gyulassy-Levai-Vitev (GLV) reaction operator approach to 
multiple elastic scatterings |8| to take into account the longitudinal momentum reduction for 
a hard jet (\/jr*/p° — 0) that propagates through dense nuclear matter due to the medium 
recoil. At the leading power approximation in 0(1/P) for a fast parton of momentum p^ = 
(P, 0j_, P) elastic scattering in nuclear matter introduces transverse momentum broadening 
without changing its "+" light-cone component. By including leading power corrections due 
to the medium response, we evaluate the parton's longitudinal momentum shift, —Ap» = 
(pj_)/2P, that ensures 3D momentum conservation in an elastic collision up to corrections of 
0({p\) 2 /P 2 ). In our derivation we introduce a boost invariant formulation of the problem 
of multiple jet interactions that generalizes the Gvulassy-Wang picture of static or massive 
scattering centers |i| employed in recent studies jj, U, Q|. 

We apply the transverse momentum diffusion results in a simple analytic model to eval- 
uate the nuclear modification of the back-to-back di-hadron correlation function C ( A<p) = 
N h ^(A(f))/N^ t M in d + Au reactions at RHIC using cold nuclear matter transport coeffi- 
cients p extracted from low energy p + A data ^(|. We find a small, weakly dependent on 
centrality, increase of the far-side (A0 > it/ 2) width o Far of C(A(f)) d A u - Such behavior is 
distinctly different from the reported disappearance of the back-to-back correlations in cen- 
tral Au + Au reactions which can be understood in terms of strong final state radiative 
energy loss 12 1 and subsequent redistribution of the lost energy in the parton system 
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II. MULTIPLE ELASTIC SCATTERING IN THE GLV REACTION OPERATOR 
APPROACH 



The GLV reaction operator formalism was developed for calculating the induced radiative 
energy loss of hard quarks or gluons when they pass through a dense medium In this ap- 
proach the multiparton dynamics is described by a series expansion in \ = J a e i(z)p(z)dz = 
L/X, the mean number of interactions that a fast projectile undergoes along its trajectory. 
Each interaction is represented by a reaction operator that summarizes the unitarized basic 
scattering between the propagating system and the medium. The summation to all orders in 



4, 



S(p+ - V2P)S 2 (p ± ) . (2) 



X is achieved by a recursion of the reaction operators and is given in a closed form in 

In this Letter we focus on the case of multiple elastic (no-radiation) scattering of a jet of 
momentum p in nuclear matter. Let M = ie ip ' z °j(p) ld R xd R be the amplitude of the parton 
in color representation R of dimension (Ir prepared at a position Zq = [zq, z^ , zq ±], where 
the light-cone coordinates are defined as z ± = (z° ± z 3 )/\^2. The unperturbed inclusive 
distribution of jets in the wave packet is given by 

For the instructive case of a normalized forward monochromatic beam (po = P\\ = P) 

d 3 N l 
dp + d 2 p± 

In Eq.(j2J) the constraint on the parton's "— " light-cone component from the p 2 = on-shell 
condition is also shown. 

In the presence of nuclear matter the multiple jet interactions are modeled by scatter- 
ing in the presence of an external non-Abelian field V J '' c (q) which satisfies the condition: 
q^V^' c {q) = 0, c being the color index. When the parton energy is much larger than the 
typical momentum scale of the medium, P 2 ^> \q 2 \, we have 

V^{q) = n» 2n5(q + ) V c (q) , g s V c (q) = v(q) T C (T) , (3) 

where g s is the strong coupling constant, the four-vector n M = 5^^ = [0,1, 0j_] and q^ = 
[q + = 0,g~,q_i_] is the momentum exchange with the medium. The phase in Eq.fjHJ) keeps 
track of the position of the scatterer relative to a fixed space-time point, e.g. the hard 
collision vertex. The color matrix T C (T) in Eq.@ represents the non-Abelian charge of the 



Mj (P) Mf(p+q) 




FIG. 1: In- medium interaction of fast almost on-shell p 2 ~ quarks and gluons with an external 
color field V^(q) T c located at position z. The momentum flow for the elastic scattering subprocess 
is shown in the diagrams. 



field radiated by a quark (or a gluon) with T = fundamental (or adjoint) representation of 
SU C (N) of dimension dj-. The Fourier transform of the non-Abelian field is given by v(q) 
and, similarly to the Gyulassy-Wang model we employ the color-screened Yukawa type 
but with Lorentz boost invariance: 



v(q) = 



Ana. 



Ana. 



(4) 



—q 2 + /i 2 q_l + /i 2 

where we have used the q + = choice of frame. This specific form of v(q) = v(q±) is 
particularly useful since in-medium interactions in both hot and cold nuclear matter are of 
finite range r^. = ^T 1 and we shall assume that Xfi ^> 1, where A is the parton's mean free 
path. 

Figure Q represents the simplest case of a quark or gluon jet with a large forward mo- 
mentum p + scattering in the medium. For initial state interactions, the elastic scattering 
amplitude for the subprocess with momentum flow illustrated in Fig. HJa) is given by: 

ij • (p + q) 



(2tt) 4 2 

dq~ d 2 q_|_ 
2tt (2tt) 2 



_(p + l) 2 + ^ 
2p+v(q ± ) 



-ig s rT c {F)) 



Jq (z-z ) 



+ e -*q±-(z±-zo±) T c (p\ T C (T) 



(5) 



where T C (F) is the color matrix in the quark representation. Before giving the corresponding 
matrix element for the gluon case in Fig. E^b) we specify the use of the light-cone gauge 
A^n^ 1 = A + = 0, a "physical" gauge for a system moving very fast along the "+" direction. 
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In this gauge we evaluate the gluon scattering amplitude: 
d A q 



M 9 el 



x 



4 [(-g s facb)(9a^(p ~ q)p + gtf(2q - p) a + 9pa(~2p ~ q) p] 
pp, | n /3 (p + q) 13 + -nP (p + qf 



(2tt 

-<7 W + 

dq~ c/ 2 q_L 
2tt (2tt) 2 

x (M 2jQ Mf) 



(p + g) 2 + «e 



(p + 5) • n 



V*(q)M 2 ,p,(p+q)M l!a (p) 



2p + v(q_ L ) 
(p + q) 2 + ie 



e iq (2-2 ) + e -iqx-(z±-zox) T C (J^\ T C (T) 



(6) 



In the derivation of Eq.(jHJ) we have used the approximations |q_J <C p + and (p±q) fl M fl (p) m 
0. The long range color interference has been neglected, allowing to factor out M" M 2)Q ,, and 
—ifacb — \T c {A)\b a is the color matrix in the adjoint representation. We note that for small 
angle scattering of high-energy partons the amplitudes M^,M^, Eqs. fl5l6|) . differ only by a 
color factor. 

The remaining q~ integral can be performed by closing the contour in the lower half- 
plane (z + — Zq < 0) and picking the contribution at q~ = q 2 L /2p + — ie, corresponding to the 
pole in the propagator of the scattered parton of momentum p + q. The factorizable elastic 
scattering amplitude becomes: 



-iv\z Q — z 



(27T) 



(7) 

Similarly, for final-state elastic scattering with the "observed" parton of momentum p M , we 
get the same scattering amplitude except the pole of the q~ integration is in the upper 
half-plane, q~ = — q 2 L /2p + + ie, and the argument of the ^-function in Eq.(jZJ) is reversed to 

In computing the elastic scattering cross section in matter one takes the average over 
the distribution of scattering centers in the transverse plane of area At, 1/A± J d 2 z± (■ ■ ■) 
which diagonalizes the squared amplitudes in the q_i_,q' ± variables J, 

T(z 



*(qx-qxM z x-zox) . ^ 

A N 



(2vr) V(q ± - q' ± ; 



(8) 



in the V 'A± ^> limit. T(zo_l) = / dz p(z,zo±) is the Glauber thickness function at 
impact parameter zo_i_ and 1/A± = T(zq±)/N with JV being the number of scattering centers. 
The corresponding differential elastic scattering cross section per unit partonic scattering 
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(T(z 0± )/N = 1) is found to be: 

da d (R,T) C 2 (R)C 2 (T) Kq ± )| 2 C 2 (R)C 2 (T) 4a 2 



(9) 



d 2 q± d A (2tt)2 d A (ql+/i 2 ) 2 

for a fast parton in color representation R and the non-Abelian field radiated from a parton 
in color representation T. For example, the color factor C 2 (R)C 2 (T) / d A = 2/9, 1/2, and 9/8 
for the scattering of qq, qg, and gg, respectively. A real nuclear medium can be a mixture 
of soft quarks and gluons with a "mean" color factor (C 2 (R)C 2 (T) / d A ) ■ For \i — 0, Eq.(jHJ) 
is effectively the small angle, t = — qfi_ and u ~ s, limit of the exact qq — > qq, qg — > qg and 
gg — > gg elastic parton-parton scattering cross sections. 

In deriving Eq.Q the 0(l/p + ) correction, oc (q 2 L /2p + — c^ 2 ± /2p + ), from the phase factors 
in Eq.flJJ) vanishes due to Eq.flHJ), <5 2 (q± — q'j_). Thus, interference effects for an energetic 
jet are largely suppressed, but they become 

case of multiple elastic scatterings the proper inclusion of unitarity corrections and their 
derivation has been discussed in Ref. [sj]. We present here the final result, noting that there 
is once again cancellation of the interference effects controlled by a phase cx (q_i_ + q'j_) 2 /2p + 
due to the forward weight, £ 2 (q_i_ + q'j_), of the virtual scattering: 



do u . c {R,T) 
d 2 q± 



^ C 2 {R)C 2 (T) |«(£j 



S 2 (q±) = -a el 6 2 (q ± ) , (10) 



d A (2tt) 

i.e. the total scattering cross section in the forward q± = direction. 

The full solution for the momentum distribution of a jet that has traversed nuclear matter 
of longitudinal extent L and opacity x = L/A = pLa e i = T(zo±)a e i can be obtained using 
the Gyulassy-Levai-Vitev reaction operator approach [8J: 



d 3 Nf(p+,p ± ). ^ X n 



dp + d 2 T>\ p~=t j j- ^— ' n! , 



1 da el (R,T) 
a e i d 2 qi± 



-qix-v 



x T-TnZ 1 p 2 , I 11 ) 



dp + d 2 p± 



2pi 



for the leading power in q^/P 2 . Eq. ()ll|) shows that in the asymptotic limit, q^/P 2 ~ 
0, elastic multiple parton scattering in nuclear matter introduces transverse momentum 
broadening to the jet without changing its "+" momentum. The average "— " momentum 
component becomes p~ = (p 2 L )/2p + from the on-shell constraint but 3D momentum is not 
conserved: |p 2 | = ( Pll ) 2 + (p 2 ) » (P 2 - (p 2 )/2) + (pi) ^ P 2 . 
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The constant "+" momentum component is a direct consequence of the leading power 
approximation, which results in the momentum constraint 5(p' + — p + ) = S(q + ) in Eq.(J3j). 
We here focus on the q^/P 2 power correction due to the medium recoil because it gives 
a leading contribution to the change of the fast parton's "+" momentum. In a simple 
model of elastic jet scattering off a non-Abelian field radiated by a parton with momentum 
k v = k u ~ = (V2P)n u the correction from the medium response replaces the momentum 
constraint at the leading power by 5(p' + — p + (l — k' + /p + )) for recoil momentum k' + . We 
find k' + = (q^/Qo + k) 2 )p + = ~(q±/y/2P) if the scattering puts the recoil parton on-shell 
(k 12 = 0). By including this longitudinal momentum reduction Eq. ljll)) is modified to: 



d?Nf(p+,p ± ) 
dp + d 2 p± 



2^ 



n=0 



i=i 



1 da el (R,T) 



Pel 



d 2 q { 



■v P1 e m<$jv2P)d p+ 



X 



d 3 N i (p + ,p ± ) 



2^ 



dp + d 2 p± 

For any initial jet flux, Eq.lfT]). the opacity series in Eq. (|12|) is most easily resummed in 
the impact parameter space (b~,b±) conjugate to (p + ,p±). We substitute in Eq.(JT2J) the 
illustrative example of an initial parton momentum distribution, Eq.(j2J), which gives 



121 



d 3 Nf 
dp + d 2 pj 



(r,bj 



: V2Pb~ 



2tt 



cxp 



X / d 2 q± 



1 da, 



a e i d 2 q A 



el -iqx-bx -ii(qiA/2P) b~ 



(13) 



While it is difficult to find a closed form for the integral in Eq. (jl3j) even for simple forms 
of the differential elastic scattering cross section, the first correction for large longitudinal 
momentum can be evaluated by expanding e~*^ q i/ v/2 " p )'' = 1 — i^(qfi_/ \/2P)b~ + • • • . The 
azimuthal integral in Eq.(|lH|l yields Jo(b±q±) and the first two terms in the expansion of the 
exponent integrated over the Yukawa potential give: 

2/i 2 



.1 b~ 



—i 



2V2PJo 



dq±q± 



dq±q 3 ± 



(fi 2 + ql) 2 
2/i 2 



^ 2 + ql) 2 



Jo( b ±Q±) = b±/iK 1 (b ± n) 
1 b- \ 



Jo( b ±Q±] 



—i 



2 v^py 



fi 2 (2K (b ± fi) -b x nK x (b x n)) 



(14) 



The key to evaluating the average p_i_-broadening and the shift in the "+" momen- 
tum of the partons is the small b±fi expansion in Eq.(fTl|): b±n Ki(b±n) = 1 — 
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(blfi 2 /2)[\n(2e-^/(b ±f i)) + l/2} + 0(bifi 4 ), K (b ±f i) = ln(2 e -^/(6 ± /i)) + C(6i/i 2 ). Keeping 
terms oc x/-* 2 ^ where £ = ln(2e~ 7J3 /(6j_/i)) > 0{1) we find: 



d 3 Nf 
dp + d 2 p^ 



1 i,„.2 C k2 +i(V2P-i2x^l 



(6-,bJ = -i-e-^ b ie T V v . (15) 

[ATI) 6 



Because of the power behavior of b 2 ± , we treat £ as approximately constant when Fourier 
transforming Eq.(15) back to momentum space: 



d 3 N f (p + , p ± ) 1 e 2 xM 2 e 

■ o 



(16) 



dp+d 2 p ± P~=^r 2vr XA* 2 £ 

The broadening of the parton beam, approximated in Eq. (jl6|) by a Gaussian, induces a 
negative light-cone component via the poles of the projectile propagators (see e.g. Eqs. (|5l6j) ) 
that ensure on-shellness at all intermediate stages. The power corrections from the target 
recoil lead to a reduction of the large "+" momentum component: p + —>■ p + — (p 2 1 )/2p + , 
where (p 2 ) = J d 2 p± p 2 ± (d 2 N f /d 2 p±) / J d 2 p± (d 2 N^/d 2 p±) = 2%/i 2 £ is evaluated from 
Eq. (fT6*j) . We note that the recoil power corrections to p~ are (9((p 2 L ) 2 /(p + ) 2 ) and hereby 
neglected. To the order to which we computed the final state distribution (p 2 ) = 2x/z 2 £, 
(p + ) = \p2P — \ 2 ^p and (p~) = 1 ■ We emphasize that proper inclusion of the target 
response ensures p° = P, \p 2 \ = (p\\) 2 + (Pj_) = P 2 energy and 3D momentum conservation 
in our final result Eq. fTo^ . 

One important consequence of the formalism that we have presented is the longitudinal 
momentum backward shift: 

dp\\ Ap,| /i 2 (20 1 ( ti 2 



— ' ( 17 ) 

dz L X^ g 2p|| \\,gj e // 2p|| 

that couples to the transverse momentum broadening, and may mimic small elastic energy 
loss if the full structure of d 3 N? / dp + d 2 p± is not observed. We note that the results on elastic 
transverse momentum diffusion presented in this Section are rather general, for example, the 
same expression is obtained for electron scattering via photon exchange or the scattering of 
a fast nucleon in nuclear matter via pion exchange. 
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III. APPLICATION TO THE BROADENING OF THE BACK-TO-BACK DI- 
HADRON CORRELATION FUNCTION 



As an application of the multiple initial and final state elastic scattering formalism elab- 
orated here we consider the nuclear induced broadening of the back-to-back jet correlations 
associated with hard QCD ah — ► cd partonic subprocesses. We will limit the discussion 
to the Gaussian 2D random walk approximation, Eq. (jl6|) . to make use of its additive dis- 
persion property. Beyond this approximation the power law corrections to the tails of the 
distribution can be numerically evaluated as in 

For p + p collisions discussion of the jet acoplanarity resulting from vacuum radiation 



is given m 



151 ]. Experimental test of the theoretical estimates is performed through mea- 



surements of the near-side and away-side di-hadron correlations in a plane perpendicular to 



the collision axis 



in. 



161 ] . In the absence of soft gluon bremsstrahlung and deviations from 



the double collinear pQCD approximation the back-to-back leading hadrons are expected 
to be coplanar. Two major effects contribute to the observed acoplanarity in p + p: the 
non-perturbative fragmentation process, where final state hadrons pick up momentum 
perpendicular to the thrust axis, and the k T dependence in the unintegrated parton distri- 
butions <f) a /h(x, ky) |l7|. Phenomenologically extracted (ky) pp vacuum broadening can be 



implemented in the pQCD hadron production formalism as in jl5j, |18( via is a normalized 
2-dimensional Gaussian. The dependence of D h / a (z,} T ) has not been discussed in the 
context of single inclusive particle production. 

Measurements of intra-jet correlations find an approximately Gaussian jet cone shape. 
If one defines { | ) to be the average particle transverse momentum relative to the 
hard scattered parent parton in the plane normal to the collision axis, it can be re- 
lated to the width a Near of the near-side (A<p < 7r/2) di-hadron correlation function 
C(A<f>) = N hl>h2 (A(j))/Nfg t ' h2 as follows: (|j Ty |) = (|p T |) sm(a Near /^). Results from the 
CCOR collaboration (consistent with preliminary PHENIX measurements) give a value 
for — 400 MeV that is roughly independent of y/s and the px of the trigger particle. 

This allows for the approximate separation of short and long distance dynamics in QCD in 
accord with the factorization theorem [19] for inclusive hadron production with momentum 
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exchange much larger than (UtJ) and { | k-j- y |) : 

E-p- oc ^2<p(x a ,k Tb ) ®<j)(x b ,k Tb ) <g> — — — ®D B / h (z e ,j T ) , (18) 

ab,c 

where "®" denotes standard integral convolution over the parton momentum fractions. 
Moreover, a lack of statistically significant broadening of the near-side di-hadron correlations 
in Au + Au can be interpreted as a signal for fragmentation sufficiently outside of the 
interaction region where the hot and dense quark-gluon plasma ^| is expected to be formed. 
The reported consistency in the s 



shape of the near-side peak of C(A0) in p + p, d + Au, and 
16] puts strong constraints on non- fragmentation and short 

n n 

formation time hadronic scattering models |2QJ that, similarly to hydrodynamics 21\, are 



Au + Au for p T > 2 GeV 
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free of jet-like correlations. 



The out of trigger plane momentum component of the far-side correlated hadron is given 
by Q: 

{\VTout\? = (\ky\) 2 + 4 i(\iTy\) 2 + 2 (|k T ,|> 2 ) . (19) 

In Eq. (|19|l xe = —Pr ' Pt trig / \pt tri g \ 2 ~ — cos(|A0|), where A0 is the angle between 



the approximately back-to-back hadrons hi,h% and (|A0|) = a/2/7t <Jp ar . Expressing also 
i\PTout\) = |Prl sin(|A0|) and eliminating it in Eq.(JTHl) for (|A0|), |p^| ps \p T tri g \ = (|Pt|) 
the following approximate relation is found: 



(\ k Ty\) = (|Pt|)cOS 



ear 



l t^[J-a Far ) -tan 2 f^) . (20) 



/ — \ \ r, UC<JJ " L \ '-'far I ^« iJ - 1 , — 

V n J \ 2 y V vr J Vv 71 " 



I , as 

3- 



In the presence of nuclear matter initial state (IS) and final state (FS) transverse momen- 
tum diffusion add a large (Aky) oc A 1 / 3 term to the vacuum parton broadening 0, Isj 
confirmed by Fermilab experiments on the nuclear A-dependence of di-jet acoplanarity 
For d + Au reactions at RHIC in the IS before the hard collision only the partons from the 
incoming deuteron scatter multiply on the nucleus: (Ak 2 ^)/^ = (fi 2 /\) e ff(L) Is . The gluon 
scattering dominated transport coefficient u 2 /X ~ 0.14 GeV 2 /fm is constrained [5] 

from ex- 
isting low energy p + A data [l(| and consistent with u 2 /X q = 0.047 ±0.035 GeV 2 /fm found 
in |23j when one takes into account that X q /\ g = Ca/Cr = 2.25. After the hard collision 
vertex in the FS both outgoing jets scatter in the medium to acquire (Ak|,) ps- Projection on 
the plane of measurement with which k-r fs forms an angle cufs, (cos 2 a fs) = 1/2, reduces 
the FS effect by a factor of two. Comparison of nuclear broadening in Drell-Yan and di-jet 
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FIG. 2: Predicted broadening of the back-to-back jet correlations in minimum bias d+Au reactions 
relative to the vacuum case is illustrated for (Ik^l) (top panel) and op ar (bottom panel). The 
strength of FS interactions is taken to be Kps = 1)2 times bigger than in the IS and represented 
by a band. Kps = 2 is a good estimate of the different amount of k^-diffusion in central versus 
minimum bias events. Preliminary p + p data is from PHENIX and the dot-dashed line estimates 
the vacuum broadening. 



production indicate that the strength of FS scattering may be bigger in the IS 22|. This 



can be modeled by Kps = (^r)e// Fs/i^eff is- Lacking precise experimental data, we nat- 
urally choose Kfs = 1 to be the default value of this parameter. The total vacuum+nuclear 
induced broadening in the plane perpendicular to the collision axis in d + Au is given by: 



(14) (14) ,.,„. V l je t ( (L) IS + 2 jets (I) K FS (£) (L) FS 

A J eff V Z / projection \ A J e// 



(21) 



To relate the effective opacity (or mean number of scatterings) \ = n = (L)/X e ff for the 
IS and the FS it is useful to look in the rest frame of the nucleus where the back-to-back 
midrapidity hard scattered partons move in the forward direction at an angle 9 ~ 7 -1 relative 
to the incident beam axis, 7 = E N /m N = y/s/(2rriN) being the Lorentz boost factor. For 
RHIC 7 = 106 and uniform distribution of the hard scatter in the nucleus xis — Xfs- 

The estimate for (|krj,|) mc = 0.75 GeV per parton is taken from (shown by dot- 
dashed line in Fig. |2J) and the di-hadron correlation function is approximated by near-side 
and far-side Gaussians. Such simplification does not take into account the large angle 
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— Preliminary PHENIX p+p vacuum broadening 

— Min. bias d+Au vacuum+nuclear (K FS =1) broadening 

— Central d+Au vacuum+nuclear (K FS = 1 ) broadening 

» Scaled STAR p+p vacuum broadening for 0.5(h + +h ) 



p+p versus minimum bias and central d+Au correlations 
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FIG. 3: Elastic broadening of the back-to-back di-hadron correlation function in minimum bias 
and central d + Au reactions relative to the vacuum p + p result (here represented by a shaded area 
to guide the eye) for (|pr|) = 4 GeV. C(Acp) is modeled via near-side and far-side Gaussians with 
away-side width evaluated from Eas. 1)20121(1 and Kps = 1- Scaled STAR correlation data for 
\{h + + h~) in p + p is also shown. 



multiple scattering and fragmentation and the growth of the background in the intermediate 
A(j) ~ 7r/2 region. The combined vacuum+nuclear matter induced IS and FS transverse 
momentum diffusion for minimum bias d + Au reactions is evaluated from Eq. (}2l"| . Noting 
that for two Cartesian components { | y | ) = a/ (k|,)/7r, in the top panel of Fig. El ~ 30% 
enhancement of ( | y | ) that measures the increased acoplanarity due to /^-broadening is 
predicted. The lower panel illustrates the increased width, o> ar , of the away-side correlation 
function C(A<f))dAu obtained solution to Ea.fZUl). 

Figure El shows the change in the shape of the far-side di-hadron correlations from trans- 
verse momentum broadening and the coupled 20% reduction of their maximum strength 
(at A<f) = 7r). In central d + Au reactions the effect of nuclear matter induced acoplanarity is 
slightly larger. We emphasize that even a significant increase in the opacity of the medium 
will lead to only a subtle change in op ar since its observable growth is no faster than ^J\- 
Triggering on very high-p^ hadrons will bias the measurement toward quark jets and may 
lead to a reduction of the estimated broadening. 

For Au+Au collisions one has to account for the IS broadening of both partons and replace 
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the FS diffusion in cold nuclear matter with quark-gluon plasma induced broadening, which 



for the Bjorken expansion scenario 



21 



(Ak 2 T ) FS = C 2 (R) 



3ira 2 s 1 dN<> 
2 A ± dy 



2£ In 



at midrapidity has the form: 
(L T ) 



T 



C2(R) = 4/3 (3) for quarks (gluons) 



(22) 



In Eq. fl22|) dN 9 /dy is the initial gluon rapidity density, {Lt) is the transverse size of the 
medium and r is the initial equilibration time 21 1. 



IV. CONCLUSIONS 

We have extended the GLV reaction operator formalism to multiple elastic scatterings 
to take into account the energy-momentum conservation for the incident parton flux and the 
recoil of the target partons along the projectile path. The resulting coupling between the 
forward and backward light-cone and transverse momentum components of the distribution 
of jets that have traversed dense nuclear matter is non-negligible for small jet energies. We 
have introduced a boost invariant formulation of the problem of multipleparton interac- 
tions that generalizes the Gyulassy-Wang model of dense nuclear matter |2| . Relaxing the 
assumptions of static or very massive scattering centers in the target simplifies the calcula- 
tion, the transition between different reference frames and also allows to analytically include 
corrections, accurate to 0(l/p + ) in Eq.(|16|). to the leading power approximation. 

In a practical application of the transverse momentum diffusion results we evaluate the 
nuclear- induced acoplanarity of the hard-scattered parton pair, Eq.(|21|). and relate it to the 
expected broadening of the back-to-back di-hadron correlation function in d + Au reactions 
at RHIC. We find that scattering in cold nuclear matter, which results in a predicted small 
20-30% Cronin enhancement also leads to ~ 25 — 30% increase of the far-side width 
op ar of C(A(t>)dAu m minimum bias d + Au reactions for { |pt| ) = 3 — 4 GeV. The centrality 
dependence of op aT is shown to be weak with ~ 35 — 40% growth of the far-side width in 
central d + Au relative to p + p. In summary, di-hadron correlations for scattering in cold 
nuclear matter are predicted in Fig. El to be qualitatively similar to the p + p case. 

Note added in proof: During the completion of this manuscript data on the nuclear 
modifications of hadron production and back-to-back di-hadron correlations in d + Au 



at \^s NN = 200 GeV became available 2J|. The moderate-px enhancement observed in 



minimum bias d + Au versus the jet quenching established in central and semi-central 
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Au + Au HQ are in quantitative agreement with the predictions in (5(. C(A<f)) d A u and 
C(A(f)) pp have the same qualitative behavior in both the near-side and far-side regions. 
Further anal ysis of the experimental data is needed to clarify the centrality and rapidity de- 
pendence P.!^ of the Cronin effect at RHIC and to quantify the broadening of the far-side 
width, a F ar , discussed here. 
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